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SUMMARY 


An  analytical  noriel  of  the  charginq  of  a  spacecraft  emitting 
electron  and  ion  beans  has  been  applied  to  the  case  of  large  spacecraft. 
In  this  ncdel ,  ionization  occurs  in  the  sheath  due  to  the  return  current. 
Charge  neutralization  of  spherical  space  charge  flow  is  examined  by 
solwinq  analytical  equations  numerically.  Parametric  studies  of  potential 
o*  larqe  spacecraft  are  performed.  As  in  the  case  of  snail  spacecraft, 
tne  inns  created  in  the  sheath  by  the  returninq  current  play  a  larqe 
role  in  deteminino  spacecraft  potential. 


INTRODUCTION 


The  potential  difference  created  between  spacecraft  qround  and 
the  annient  olasna  durinq  the  ejection  of  a  bean  of  electrons  fron  a 
soijn-jino  rocket  payload  in  the  ionosphere  (Pef.l)  has  been  found  to  be 
<‘"Ch  less  than  had  oriqinall.y  been  theoretically  predicted  (Ref.?). 

To  determine  the  reasons  for  this  limited  potential  difference,  large 
vacuum  cha  ber  tests  were  conducted  in  which  electron  and  ion  currents 
were  elected  from  a  p'ayl'oad  into  a  simulated  ionosphere. 

As  a  plausible  explanation  to  the  observed  current  voltage  beha¬ 
vior,  sheath  ionization  models  (Refs.  3,4)  for  snail  spacecrafts  have 
beer-  studied.  When  an  electron  beam  is  emitted  from  a  spacecraft,  ambient 
electrons  are  attracted  by  the  charqed  spacecraft  (Ref. 5).  They  collide 
with  the  neutral  atmospheric  molecules  in  their  paths,  and  may  be  ener¬ 
getic  enouuh  to  ionize  the  neutrals  to  fo-m  new  electrons  and  ions  (RefL 
6).  These  newly  created  charqes  alter  the  space  charge  current  arriving 
at  the  spacecraft  shiftinq  the  potential  to  a  lower  value.  The  bean 
electrons  are  assumed  to  he  enerqetic  enough  to  leave  the  spacecraft 
completely  an*  to  play  a  neqliqible  role  in  the  Ionization.  This  mechanism 
is  capable  of  explaining  the  non-nonotonic  current-voltage  behavior  ob¬ 
served. 
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In  this  paper,  we  apply  the  sheath  Ionization  model  to  the  case 
of  large  spacecrafts  In  the  Ionosphere.  In  particular,  it  Is  Important 
to  find  out  whether  the  non-monotonic  current-voltage  behavior  during 
electron  beam  emissions  would  still  be  present  for  large  spacecrafts. 
Details  of  the  method  are  given  and  followed  by  a  discussions  of  results. 


SYMBOLS 

electron  charge 
electric  field 
beam  current 
mass  of  electron 
mass  of  ion 

density  of  ambient  electrons 
density  of  Ionization  Ions 
density  of  Ionization  electrons 
probability  of  ionization 
radius  of  spacecraft 

radial  position  measured  from  the  center  of  spacecraft 
radial  position  used  as  Integration  variable 
radius  of  sheath  measured  from  the  center  of  spacecraft 
velocity  of  an  electron  In  the  sheath 


thermal  velocity  of  an  ambient  electron 
sweep  velocity 

mean  free  path  of  electron  neutral  collision 
electric  potential 
space  charge  density 
permittivity  of  apace 
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MATHEMATCAL  FORMULATION 


The  method  of  approach  used  Is  to  study  an  analytical  'plasma  probe* 
model  (Refs.  7,8),  with  space  charge  flow  of  electrons  accelerating  through 
the  sheath  suroundlng  a  spherical  ' probe' ,  which  represents  a  spacecraft  In 
an  lonlzable  plasma  environment.  Magnetic  field  effect  is  ignored  In  this 
model. 


Hie  beam  Is  assumed  to  be  energetic  enough  to  leave  the  spacecraft 
completely,  and  is  not  stopped  by  Its  own  space  charge  at  all.  As  the 
beam  electrons  leave,  the  spacecraft  becomes  charged  oppositely.  A  pola¬ 
rization  region  (sheath)  Is  formed  In  the  vicinity  of  the  spacecraft.  In 
our  model,  ions  are  assumed  depleted  due  to  charge  repulsion  Inside  the 
sheath  (Figure  1). 


The  depletion  radius  r0  will  be  defined  by  the  balance  of  the  out¬ 
going  beam  current  with  the  incoming  ambient  current.  For  a  beam  current 
If,,  the  depletion  radius  r0  la  determined  by 


V4'  ro2°eevth 


(1) 


where  vth  is  the  thermal  velocity,  and  la  the  number  density  of  am¬ 
bient  electrons.  Some  typical  values  of  sheath  radius  as  calculated  by 
means  of  eq(l)  are  shown  in  figure  2. 


The  potential  +  at  any  point  Inside  the  sheath  is  governed  by 
Poisson's  equation  : 


V  2  4 


(2) 


w*>era  p  is  the  apace  charge  density,  and  c0  Is  the  permittivity  of 
empty  space. 
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Spherical  Symmetric  System 


To  simplify  the  geometry,  we  assume  spherical  symmetry  In  the  space¬ 
craft  and  sheath  system.  Equation  (2)  becomes  simply  a  radial  equation  : 
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r 
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♦  (r)  \ 
3  r  / 


P  (r) 


(3) 


where  the  gradient  of  the  potential  t  gives  the  electric  field  E  : 


3*  (r) 

-  -  -  E(r)  (4) 

3  r 


Taking  into  account  the  electron  and  Ion  pairs  created  as  a  result  of  Ion¬ 
ization,  the  charge  density  p  at  say  point  r  In  the  sheath  Is  given  by  the 
sum  of  charge  densities  (figure  3): 


p  (r)  -  e  [n+(r>  -  n"(r)  -  ne(r)] 


(5) 


where  ne  is  the  return  current  (primary)  electron  density,  while  n+  and  n" 
are  the  ionization  ion  and  electron  densities  respectively,  due  to  return 
current  electron  collisions  with  neutrals. 


The  Ionization  electron  density  n-(r)  Is  due  to  all  Ionizations 
that  occur  outwards  of  r,  and  the  density  n+(r)  of  Ions  at  r  Is  due  to  all 
Ionizations  that  occur  Inwards  of  r.  Thus, 


n”(r) 
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r2 
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and 


n+(r) 
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where 


f  dn  1 


>  ”1P(ve(r’)]ne(r')ve(r<) 


(8) 


Numerical  Method 


To  solve  the  system  of  equations  (3)  to  (8),  one  divides  the  space 
of  the  sheath  into  N  concentric  shells,  and  sets  up  N  equations  for  the  N 
unknowns  <J>  i  (see  figure  4).  In  view  of  the  complexity  of  the  ionization 
terms  in  equations  (6)  and  (7),  it  is  Impossible  to  solve  these  equations 
exactly.  Instead,  one  seeks  the  approximate  solutions  that  minimize  a 
function  F,  the  mean  square  of  f^,  constructed  from  the  radial  Poisson 
equation  (eq.3)  for  the  i-th  cell,  where  i*l,  ...»  N. 


*  (r  E)j^j_(r  E)j*  — —  (r  p  (^,  t  (9) 

Eft 


where  the  electric  field  E  (eq.4)  is  constructed  in  a  finite  difference 
scheme  : 


♦  i  •  ♦  1+1  "  A  rtE^Ei-H+Ej+j)^*  (10) 
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The  numerical  method  used  to  solve  equations  (8)  to  (10)  Is  the 
standard  Newton- Ralphson  laethod  of  Iteration  : 


=  E,U> 


f1(E1(J>,...EMU>) 
a  f1(E1^)....EM^,)/5  El 


(11) 


A  set  of  trial  solutions  Is  used  to  start  In  the  Newton-Ral  phson 
Iteration  process,  and  a  convergent  set  of  solutions  Is  seeked,  for  each 
set  of  Input  parameters  such  as  beam  current,  ambient  electron  density, 
ancient  electron  temperature,  mean  free  path,  and  spacecraft  radius. 


RESULTS  AND  DISCUSSION 


Figure  5  shows  the  computed  results  of  spacecraft  potential  as 
a  function  of  electron  beam  current  for  various  electron  density,  elec¬ 
tron  temperature,  and  mean  free  path.  The  non-monotonlc  behavior  of 
potential  current  curves  shows  up.  At  low  currents,  the  potential  In¬ 
creases  with  beam  current.  When  the  current  Increases  further.  Ioniza¬ 
tion  occurs  Inside  the  sheath.  The  potential  then  turns  around  as  the 
current  of  the  electron  beam  increases. 


The  ton  and  electron  charges  created  by  ionization  alter  the 
behavior  of  the  space  charge  flow,  originally  governed  by  the  single 
charged  Poisson  equation.  The  potential  turns  to  a  lower  value,  and 
stays  approximately  constant  as  current  further  increases. 


At  this  stage,  the  potential  profile  as  a  function  of  radial 
distance  show  locally  flat  gradient.  This  Is  due  to  ions  created  In¬ 
side  the  sheath  not  being  able  to  move  out  quickly  due  to  their  heavy 
masses.  If  a  local  Ion  charge  build-up  forms  a  potential  hump,  ion 
motion  would  be  twoways,  and,  the  theory  would  then  break  down. 


To  overcome  this  difficulty,  a  sweep  velocity  vs  Is  added  to 
the  ions.  Eq. ( 7 )  becomes 


n*(r) 


|Jjn^  r'2  dr' 

[2elf(r)-*(r' )!/■,*  +  v$2]1/2 


(12) 
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It  is  argued  that  the  motion  of  a  spacecraft  relative  to  Its  plasma 
environment  can  provide  such  a  sweep  velocity  vs  (eq. (12)).  The  value  of 
vs  is  of  the  order  of  spacecraft  velocity  and  Is  an  arbitrary  Input  to 
the  computation.  However,  at  a  higher  current,  a  potential  hump  again 
shows  up,  and  the  computation  fails  to  converge.  The  technique  breaks 
down.  It  is  conjectured  that  twoway  space  charge  flows  should  be  accom¬ 
modated  when  a  potential  hump  appears. 


For  increasing  spacecraft  radii,  the  non-monotonic  current-voltage 
behavior  still  persists  (see  figure  6).  However,  increased  spacecraft 
radius  lowers  the  maximum  spacecraft  potential  induced  by  beam  emission. 
Also,  the  amplitude  of  the  difference  between  the  maximum  potential  and 
the  minimum  (beyond  the  turn-around)  diminishes.  Figure  7  shows  a  plot 
of  the  envelope  of  maximum  and  minimum  potentials  for  various  spacecraft 
radii. 


For  a  given  beam  current  1^  (  see  eq.(l)  ),  the  sheath  surface 
area  remains  constant  and  is  unaffected  by  the  increase  in  spacecraft  ra¬ 
dius.  The  sheath  thickness  (defined  as  the  sheath  radius  minus  the  space¬ 
craft  radius),  however,  diminishes.  As  a  result,  a  lower  spacecraft  po¬ 
tential  is  sufficient  to  attract  ambient  electrons,  through  the  sheath, 
for  the  compensation  of  electron  beam  current  leaving  the  spacecraft. 


Beyond  the  turn-around  point  in  a  current-potential  curve,  the 
minimum  potential  is  limited  by  the  minimum  energy  required.to  ionize 
a  neutral  molecule  in  the  atmosphere.  Since  such  a  minimum  energy  is 
generally  of  the  order  of  20  eV.  (Ref. 6),  the  minimum  potential  in  a 
current-potential  curve  is  expected  to  approach  about  20  eV.  aymptoti- 
cally,  depending  on  the  model  of  ionization  used.  For  the  same  reason, 
if  the  maximum  potential  induced  by  beam  emissions  is  below  about  20  eV., 
no  non-raonotonic  behavior  is  expected. 


Figure  7  shows  the  calculated  envelopes  of  the  maximum  and  mini¬ 
mum  (beyond  turn-around)  potentials  for  various  spacecraft  radii,  in  a 
given  ambient  environment.  The  amount  of  ionization  becomes  very  small 
as  the  sheath  potential  approaches  the  minimum  ionization  potential. 
The  amplitude  of  the  potential  drop  beyond  the  turn-around  also  approa¬ 
ches  the  value  of  minimum  ionization  energy. 


There  is  another  critical  beam  current,  which  manifests  itself  for 
large  spacecrafts,  but  not  for  small  ones.  This  current  is  determined  by 
equating  the  sheath  radius  to  the  spacecraft  radius.  If  the  sheath  radius 
is  too  small,  the  spacecraft  would  receive  enough  ambient  electrons  to 
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compensate  beam  emissions  without  being  charged  up.  The  potential  of 
the  spacecraft  is  that  of  natural  charging,  in  this  case.  Beyond  this 
critical  current,  the  beam  emission  is  able  to  swing  the  spacecraft  to 
an  opposite  potential  and  hence  control  the  charging  of  the  spacecraft. 
This  phenomenon  shows  up  in  the  calculations  (see  figure  6). 


In  the  model  studied,  as  the  radius  of  a  spacecraft  increases, 
three  regimes  of  physical  behavior  can  be  identified.  Figure  8  shows 
these  regimes  clearly.  The  potential  versus  spacecraft  radius  curve 
is  relatively  flat  in  the  small  radius  regime.  This  is  the  regime  in 
which  saturated  ionization  occurs,  i.e.  this  is  the  regime  beyond  the 
minimum  potential  in  a  current-voltage  curve.  The  second  regime  is  char¬ 
acterized  by  the  presence  of  the  potential  maximum,  which  is  the  main 
feature  of  non-montonic  behavior.  The  third  regime  occurs  when  the  space¬ 
craft  is  so  large  that  its  radius  exceeds  the  sheath  radius  (measured 
from  the  spacecraft  center)  for  a  given  current.  The  beam  loses  its  con¬ 
trol  of  the  spacecraft  potential,  and  natural  charging  dominates. 
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Figure  3.  Ionization  Pair-Creation  in  Sheath 


Figure  4.  [^composition  of  Sheath  into  Shells 
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igure  7.  Envelope  of  the  Potential  Extrema  in  Figure  6. 


Hqure  8.  No n-r. ont on ic  Rebavior  of  Spacecraft  Potential  as  a  Pmction  of 
Spacecraft  Radius,  for  a  Given  Electron  Ream  Current. 


